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Abstract. The hyperfine interaction between the electron spin and the nuclear spins is one of the main sources of decoherence 
for spin qubits when the nuclear spins are disordered. An ordering of the latter largely suppresses this source of decoherence. 
Here we show that such an ordering can occur through a thermodynamic phase transition in two-dimensional (2D) Kondo- 
lattice type systems. We specifically focus on nuclear spins embedded in a 2D electron gas. The nuclear spins interact with 
each other through the RKKY interaction, which is carried by the electron gas. We show that a nuclear magnetic order at finite 
temperature relies on the anomalous behavior of the 2D static electron spin susceptibility due to electron-electron interactions. 
This provides a connection between low-dimensional magnetism and non-analyticities in interacting 2D electron systems. We 
discuss the conditions for nuclear magnetism, and show that the associated Curie temperature increases with the electron- 
electron interactions and may reach up into the millikelvin regime. The further reduction of dimensionality to one dimension 
is shortly discussed. 
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INTRODUCTION 

A major source of decoherence of electron spin qubits 
is the hyperfine coupling of the electron spin with the 
surrounding disordered nuclear spins [1]. If we want 
to control and eventually eliminate this source of deco- 
herence, it is essential to fully understand the behavior 
of both the electron spin and the ensemble of nuclear 
spins. In this text we discuss such a fundamental aspect. 
We address the question whether the nuclear spins can 
achieve order through a (ferro)magnetic phase transition. 
If they do the decoherence source of the hyperfine in- 
teraction is massively reduced [2]. We focus specifically 
on GaAs-based semiconductor heterostructures confin- 
ing two-dimensional electron gases (2DEGs). Such sys- 
tems serve as the parent system for the single electron 
quantum dots defining the spin qubits |3[] . Yet we empha- 
size that the described physics remains valid for general 
Kondo-lattice systems with interacting electrons. It turns 
out that the anomalous properties of the 2DEG, result- 
ing from electron interactions, are crucial for the nuclear 
magnetic order. The exposure to follow is an overview of 
our recent work published in ]4l |5|] . 

Even though we shall be concerned mainly with fully 
translationally invariant 2DEG, let us introduce the prob- 
lem by considering a single electron spin, confined to a 
quantum dot. The important interaction discussed here is 
the hyperfine coupling between the electron spin on the 
dot, S = (S x ,S y ,S z ), and the surrounding lattice of nu- 
clear spins I, = (If,lj,If) (i is the index for the lattice 



site at position r,). The interaction Hamiltonian can be 
written in the form 



H% p = £A, S • I, = J> ;.V7 ; + S+ir + S-l? 



(1) 



with S ± = S x ± iS y and if = If ± if, and where A, w 
A|i//(r,)| 2 withA a proportionality constant and y/(r,) the 
wavefunction of the confined electron on the quantum 
dot. The number of nuclear spins is large, typically of 
order 10 5 , and so the nuclear spins generally act as a 
disordered bath on the electron spin. While If has a role 
similar to an external magnetic field, the last two "flip- 
flop" terms in Eq. (Q3 flip the electron spin and lead to 
the decoherence. 

This decoherence source can be largely suppressed 
when the nuclear spins order, ferromagnetically or dif- 
ferently, which effectively suppresses the flip-flop terms 
in the Hamiltonian [2J]. There are two ways of polarizing 
the system, which we shall call the extrinsic (or dynamic) 
and the intrinsic (or thermodynamic) polarization of the 
nuclear spin system. The extrinsic polarization consists 
in an active manipulation of the nuclear spins by the ex- 
perimentalist. Several methods to do this have been pro- 
posed and partially experimentally realized: The devel- 
opment of quantum control techniques that effectively 
lessen or even suppress the nuclear spin coupling to the 
electron spin 
distribution 
the nuclear spins fljj, 
the spin decay time by one order of magnitude through 



] ; the narrowing of the nuclear spin 
|; or the dynamical polarization of 
ImHIlBl]. Yet in order to extend 



polarization of the nuclear spins, a polarization of above 
99% is required [2], quite far from the best result reached 
to date in quantum dots, which is about 60% 

It is possible, however, that a full polarization is 
achieved intrinsically as well, i.e. through a thermody- 
namic phase transition to, for instance, a ferromagnetic 
state. This is our main topic here. In what follows we 
give a qualitative, physical account to this possibility by 
introducing step by step the model, the necessary condi- 
tions, and the results. For details we refer to Jj,|5|]. We 
first indicate how to obtain from a microscopic model 
an effective Hamiltonian for the nuclear spins only. We 
show that the conditions of the Mermin- Wagner theorem 
are not met so that long range order is not forbidden 
in the 2D system. With a simple mean field theory we 
can then see under which conditions ferromagnetic order 
is possible and allows us to identify the first important 
temperature scale Tmf- A refinement leads to the second 
important temperature scale T*, which depends on the 
electron-electron interactions in the 2DEG. The calcula- 
tion of the susceptibility for interacting electrons allows 
to estimate T* and on the stability of the magnetic order. 
It turns out that nonanalytic corrections to the electron 
spin susceptibility are crucial. We discuss the possible 
forms of these corrections, allowing us to conclude with 
some numerical estimates for a possible nuclear order. 

MODEL AND EFFECTIVE MODEL 

In order to discuss thermodynamics we have to shift our 
point of view from electrons in quantum dots to a fully 
translationally invariant 2DEG as it can be obtained, for 
instance, in a GaAs heterostructure. The Hamiltonian of 
such a system can be written as 

H = H e i+H hyp , (2) 

where H e / describes the Hamiltonian of the interacting 
electron gas, and 

i/ /iyp =£AS,-I, (3) 

i 

is the hyperfine interaction between the electron and nu- 
clear spins. We have chosen here a tight binding, Kondo- 
lattice formulation of the problem, where S; = (Sj,Sf, 5?) 
is the operator of an electron spin in a Wannier state cen- 
tered at lattice site i and the I, = (If,lJ,If) are the nuclear 
spin operators as before. For GaAs we have 1 = 3/2. No- 
tice that A now is position-independent due to the trans- 
lational symmetry. 

In Eq. © we have not included the direct dipolar inter- 
action between the nuclear spins. This interaction has the 
smallest energy scale in the system, E&& ~ 100 nK Hill . 
It is much weaker than the effective nuclear spin inter- 
action discussed below and, in particular, E^ is much 



smaller than typical experimental temperatures. This al- 
lows us to entirely neglect the direct dipolar interaction. 

For GaAs, A w 90 jxdV |3, which compares to typi- 
cal Fermi energies of Ep « 10 meV 11511 . The small ratio 
A/Ep r-j 10~ 2 implies a separation of time scales. Elec- 
tron relaxation times are much shorter than typical time 
scales of the nuclear spins. This allows us to decouple 
the systems and focus on the magnetic properties of the 
nuclear spin system alone, where the effective spin-spin 
interactions are carried through the response of the equi- 
librium electron gas to local magnetic excitations, i.e. the 
electron spin susceptibility. 

Technically, the first step is to reduce the still quasi- 
2D (due to the finite thickness of the 2DEG) problem 
to a true 2D problem. Since the electrons are confined 
in a single mode in the direction orthogonal to the 2D 
plane, the nuclear spins along a column in this direction 
are all almost identically coupled to the electrons and 
so effectively locked in a ferromagnetic alignment and 
behave like a single (effectively large) nuclear spin. The 
problem becomes, therefore, truly 2D. 

The effective Hamiltonian describing this situation can 
then be obtained, for instance, through a Schrieffer- Wolff 
transformation followed by the integration over the elec- 
tron degrees of freedom j4j, |5J . It is given by 

ija qa 

where N is the number of sites in the system, a =x,y,z 
and the lattice indices i,j run over the two-dimensional 
lattice with site positions r;,r.-. Furthermore 

J«(T i -T j ) = ~ X a (T i -Tj) (5) 

is the effective Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction [16], with n s = a~ 2 the nuclear 
spin density and 

X a (r i -r J ) = -Lj~dt{{S?(t),Sj(0)})e-*>, (6) 

is the static electron spin susceptibility (77 > is in- 
finitesimal). The Fourier transforms are defined as /" = 
£. e /r, qja and ja = J dre- iT V a (r), and N is the num- 
ber of sites in the system. We shall henceforth consider 
only isotropic electron systems, allowing us to drop the 
a index in % and J. 

Much of the magnetic properties of the nuclear spins 
depends, therefore, directly on the shape of J q , i.e. on 
the electron susceptibility #(q). Below we will first in- 
vestigate which features of #(q) are required such that 
nuclear ferromagnetism is stable. Then we shall see that 
electron interactions can indeed lead to such a behavior. 

Yet before starting we have to comment on the 
Mermin- Wagner theorem lfl7ll . which states that long 



range order in Heisenberg-like models in 2D is impos- 
sible, provided that the interactions are sufficiently short 
ranged. The RKKY interaction, however, is long ranged. 
But it is also oscillatory, and it has been conjectured 
recently lfl8ll that the Mermin- Wagner theorem extends 
to RKKY interactions that are carried by noninteracting 
electrons. Below we find a direct confirmation of this 
conjecture, indicating that electron interactions in addi- 
tion to the long range character of 7 q play the crucial role. 



MEAN FIELD THEORY 

As a first (naive) approach, we can look at the problem 
of nuclear magnetism on the mean field level, similar to 
the approach used by Frohlich and Nabarro (FN) for bulk 
metals more than 60 years ago [19]. We skip here the ex- 
plicit calculation as it is a standard Weiss mean field cal- 
culation. Instead, let us clearly state the main assumption 
behind FN's approach: If we look at the Hamiltonian 
we see that the energy can classically be minimized if the 
7™ align in a single spatial Fourier mode q correspond- 
ing to the maximum of J q . If this maximum is reached 
at q = 0, the ground state is ferromagnetic. FN implicitly 
assumed that the physics is entirely determined by this 
maximum energy scale Jo and investigated the Hamilto- 
nian HfN = —SYdj^-i ' I/> f° r which the mean field the- 
ory is exact. 

We have to retain two important points from this the- 
ory: On the one hand, the ground state depends on the 
maximum of the 7 q , and a ferromagnet is only possible 
if it is reached at q = 0. A q ^ implies a different mag- 
netic order, e.g. a helimagnet. On the other hand, this the- 
ory depends on a single energy scale, max q y q (e.g. Jq), 
which we associate with a temperature Tmf through 



x,(g)/\xs(o) 



ksTMF = max/, 



(7) 



where kg is the Boltzmann constant. Since this is the only 
energy scale in the system, all thermodynamic quantities 
must directly depend on it. For instance, the Curie tem- 
perature is given by lfl9ll 



jMF _ 



1 



-Tmf- 



(8) 



REFINEMENT 



The mean field theory is, however, inconsistent with 
the RKKY interactions for noninteracting electrons. In 
the noninteracting case the electron spin susceptibility 
is exactly known and is given by the Lindhard function 
l20ll . which is constant for < |q| < 2k F [see Fig.Q](a)], 
where kp is the Fermi momentum. There is therefore no 
well defined q at which J q is maximum. "Guessing" a 




FIGURE 1. Electron spin susceptibility %(q) for (a) nonin- 
teracting electrons in 2D (Lindhard function), (b) interacting 
electrons such that a nuclear ferromagnet would be stable. In 
case (b) the interactions introduce a new scale T* associated 
with the slope of the curve at small q [see Eq. <10H, represented 
by the thick line next to the curve. 



ground state about which we can proceed with a mean 
field theory as above is no longer possible. 

In fact, it is straightforward to show that such a ground 
state cannot be stable. For this purpose we look at fluctu- 
ations about an ordered ground state. In Heisenberg type 
systems the lowest lying excitations are magnons, collec- 
tive long ranged spin wave excitations. The calculation of 
the magnon dispersion about a ferromagnetic order is a 
standard textbook exercise (see e.g. 02111 ) and leads to 



hco„ 



2I(J -J q ). 



(9) 



We shall henceforth assume that J q is independent of 
the direction of q and write J q , %(q), as well as co q . 
The fact that J q is constant for < q < 2kp means 
that there is a continuum of magnon excitations at zero 
energy co q = [see Fig.|2](c)], and every magnon slightly 
decreases the global magnetization. The assumed ground 
state magnetization cannot be stable, the nuclear spin 
system is disordered. This is a direct illustration of the 
extension of the Mermin- Wagner theorem conjectured in 
& 

From Eq. (0 we see, however, that a ferromagnetic 
ground state becomes stable if co q [i.e. x(<l)] increases 
monotonically with q [Fig. Q] (b) and Fig. [2] (d),(e)]. 
As we shall see below, nonanalytic corrections to %{q) 
by electron-electron interactions can indeed lead to a 
linear increase in q. The main effect is illustrated in 
Fig.[TJ(b). The electron-electron interactions modify the 
RKKY interaction J q in that they introduce a new energy 
scale T*, which characterizes the shape of J q . Stronger 
electron-electron interactions lead to a larger slope of 
the linear increase of x(q) at small q. Using this slope 
and 2kp as the only inverse length scale available for the 
electron gas, the new energy scale must be set by the 
quantity 



k R T* 



{2k F 



dq 



(10) 



q=0 



The physics is, therefore, no longer dominated by the 
single scale Tmf- Thermodynamic quantities such as the 
magnetization per site m(T) or the critical temperature 
T c must be a function of these available scales: m(T) = 
m(T MF ,T*;T),T c = T C (T MF ,T*), etc. 

We can again use the magnon description to shed more 
light on this dependence. The average magnetization per 
site can be written as 

m =i-Ly ! =i-± [ ^1 

»^eH/W-l 2nJ etim q /k B T_ l > 

(11) 

with a the nuclear lattice constant. The summa- 
tion/integration runs over the first Brillouin zone of 
the nuclear system. We see that this integral converges 
if (Oq °< q for q — > 0, and so the linear corrections to the 
susceptibility are essential for the existence of a finite 
critical temperature. The magnon integral is dominated 
by this linear behavior up to T ~ T*. For these temper- 
atures, we can explicitly calculate the magnon integral 
and obtain 

m = l(l-T 2 /T 2 ), (12) 

with 

r = -^/^r~^r, ( i 3) 

Ikf \ 7t a 

where Xf — 2n/kp is the Fermi wavelength. 

We stress that the calculation leading to Eq. ( fT2l is 
valid only in the limit T < T*, where the spin waves form 
a dilute gas. At higher temperatures the number of spin 
waves increases and their wavelengths become shorter, 
leading to a breakdown of the spin wave theory. The tem- 
perature Tq in Eq. (fT2l has to be interpreted accordingly: 
We see that Tq sets a characteristic temperature scale for 
the magnetization m(T). Further corrections to Eq. (fT2l . 
and mainly the dependence on Tmf, are of exponential 
form, ~ q-Tmf/t > an( j un iq Ue iy determined by the 
spin wave modes beyond the validity of the theory. In 
addition, from Fig. [TJ we see that T* generally (i.e. for 
not too weak interactions) is comparable to Tmf, which 
is proportional to the maximum of %{q). This means that 
up to temperatures T ~ Tmf, the magnetization is essen- 
tially independent of the Tmf scale. 

If we conjecture on this basis that Tmf indeed is in- 
significant for the thermodynamics, then Tq must set the 
scale for the Curie Temperature T c , 

T C ~T ~—T*. (14) 
a 

In contrast to the mean field result ©, this estimate is 
consistent with T c — > for noninteracting electrons. The 
prefactor Xp/a^ 10 2 (in GaAs) is a consequence of cou- 
pling the electron system (with length scale Xf) to the 
nuclear spin system (with length scale a), and leads to a 
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FIGURE 2. Possible shapes for the magnon dispersion a q 
[Eq. J5J], equal to the shifted and rescaled electron spin suscep- 
tibility x{q). Noninteracting electrons lead to the curve (c), ob- 
tained from the Lindhard function. Electron interactions mod- 
ify the shape of this curve and can in principle lead to any 
of the curves (a)-(e). Ferromagnetic order is unstable in cases 
(a) and (b) where the magnon dispersion would become nega- 
tive; a helical magnetic order becomes here possible. Nuclear 
magnetism is unstable in case (c), which has a continuum of 
excitations with zero energy. A nuclear ferromagnet is stable 
for cases (d) and (e). 



strong increase of the characteristic temperature Tq com- 
pared with T* . A further discussion of the renormaliza- 
tion of Tq can be found in J^t] . 

In order to give a numerical estimate of Tq we need 
to investigate the interactions between electrons in the 
2DEG. 



INTERACTION CORRECTIONS TO THE 
SUSCEPTIBILITY 

Interactions and correlations between particles become 
increasingly important in low dimensions due to a reduc- 
tion of the available phase space for particle scattering. It 
is not too surprising, therefore, that electron-electron in- 
teractions in the 2DEG lead to deviations from the ther- 
modynamics of standard Fermi liquids. Such deviations 
have attracted some attention over the lastyears theoret- 
ically illllJMIIllIil^lMIMIMIIllQ and ex- 
perimentallv 13311. Of specific importance for our case are 
Refs. 112.61 1271 12811 . where explicitly self-energy correc- 
tions to the susceptibility were calculated for a screened 
Coulomb interaction U (x — y) = U8 (x — y) in 2D. Non- 
analytic behavior appears at second order perturbation in 
U and leads to 

Sx(q)=X(q)-X(0) = -q 4lX f}f sl2 , q^2k F , 

5%Kf 

(15) 

where F s = —Urn/An is the bare 2kp backscattering ver- 
tex and m the effective mass. We emphasize the nonana- 
lytic behavior on the modulus q = |q|, which cannot be 



derived within a standard Fermi liquid theory. This lin- 
ear |q| dependence is indeed the necessary dependence 
for a nuclear ferromagnet as discussed above. The prob- 
lem here is, however, the sign of this correction: If we 
feed this % (q) back into the spin wave spectrum (0, we 
see that co g becomes negative for g ^ 0, meaning that the 
assumption of ferromagnetic order is incorrect and that 
such an order is in fact unstable. In such a situation a dif- 
ferent, helical order can be possible, and this interesting 
possibility is discussed in |5fl- Here we focus instead on 
a further renormalization of the nonanalytic correction, 
which can reestablish the ferromagnetic order. 

We restate the comment made at the beginning of this 
section that correlation effects are important in low di- 
mensions. A result from perturbation theory, such as Eq. 
(fT5T >. undergoes further renormalization by higher order 
processes and can so change considerably its shape. This 
indeed can happen if we push further the diagrammatic 
calculation of 11261 |27t 28] and include the full summation 
of selected classes of diagrams such as, for instance, the 
Cooper channel renormalization of the two-particle scat- 
tering vertex 03411 . The Cooper channel renormalization 
for x ( T ) has been considered recently in |29[ l30l l3ll l32ll . 
The effect of this renormalization on % (q) has been es- 
timated in [51] and very recently explicitly calculated in 
lHH]. The obtained corrections to %{q) are nonuniversal 
and depend on detailed cutoff scales of the renormaliza- 
tion such as Fermi energy, temperature, level spacing, 
etc. The stability of the ferromagnetic (or helical) phase 
seems, therefore, to depend on a quantity which is diffi- 
cult to control. 

Yet we must note that this is the result of a per- 
turbative renormalization, the summation over selected 
classes of diagrams, as well as the result of a screened 
short ranged electron-electron interaction. Alternative 
approaches lead to more predictable results: If we con- 
sider long ranged Coulomb interactions within a local 
field factor approximation (which is a semi-empiric gen- 
eralization of the RPA approximation 12011 ). the spin 
wave spectrum is always positive and a ferromagnet is 
stable. Such effective theories, however, have eventu- 
ally the same difficulty of control as the summation over 
classes of diagrams. Due to this, it seems so far that the fi- 
nal determination of the shape of the susceptibility prob- 
ably has to rely on experiments and numerics. The lat- 
ter two approaches should then probe not only the slope 
dx/dq at q w but also %{q ~ 2kp), which is important 
to discriminate between the different scenarios shown in 
Fig.|2]by the curves (a) and (b) on the one hand, and by 
(e) and (f) on the other hand. We stress that numerics 
should directly target % (q) and not, as for instance done 
in ll36l l37ll . the local field factor. The relation between 
the latter and %{q) is actually singular at q = 2k,F, which 
amplifies the noise in the Monte Carlo data of j36l l37ll 
and makes a conclusion on % (q) unreliable. 



Assuming, however, that we have found a window 
in which a nuclear ferromagnet is stable, we can esti- 
mate Tq from the different calculation schemes. A de- 
tailed discussion is given in [5]. We find that remarkably 
all schemes provide comparable values. As anticipated 
much depends on the strength of the electron interac- 
tions, which can be quantified by the commonly used 
dimensionless parameter r s (see e.g. l20ll ) expressing 
roughly the ratio of Coulomb over kinetic energies of the 
electrons. In the 2DEG r s scales with the electron density 
n e and kp as r s «1/ Jn~ e « kp 1 . Values up to r s ~ 8 can 
be reached experimentally nowadays. With increasing r s , 
the scale Tq is enhanced through two main effects: First, 
kp 1 increases linearly with r s . Second, |^(0)|, which is 
essentially the Pauli susceptibility at small r s , increases 
linearly with r s [20]. We moreover note that larger r s 
drive the system closer to the ferromagnetic Stoner insta- 
bility, which would occur at r s ~ 20 [38]. At this instabil- 
ity %(0) would diverge. The proximity to the instability 
provides an additional prefactor enhancing \% (0) | . Let us 
note that the increase of |^(0)| enhances both, T* and 
Tmf- For r s ~ 5, we then obtain 7b ~ 0.3 — 0.4 mK, and 
for r s ~ 8 the larger Tq ~ 0.7 — 1 mK, where the spread 
is due to the different analytic approaches ||5J. 

ONE-DIMENSIONAL SYSTEMS 

With a further reduction of dimensionality electron- 
electron correlations become even more important. One- 
dimensional conductors of interacting electrons, such as 
quantum wires or carbon nanotubes, form a Luttinger liq- 
uid rather than a Fermi liquid. Accordingly the shape 
of the RKKY interaction J q between the nuclear spins 
in these systems changes drastically. J q is dominated 
by backscattering processes at q — 2k p, and the nuclear 
spins order in a helical phase with this wave vector. 

In contrast to the 2D case, the feedback of the ordered 
nuclear field on the electrons is now crucial. It leads to a 
spontaneous restructuring of the electron wave functions, 
i.e., to an order in the electron system as well. The feed- 
back stabilizes this order, and the critical temperature T c 
can increase by several orders of magnitude compared 
with the case where the feedback has been neglected. 

We have performed a detailed study of this effect in 
l39ll for the example of single-wall carbon nanotubes 
made from the 13 C isotope (which has a nuclear spin 
7=1 /2). Such nanotubes have become available very 
recently ll40l l4ll l42tl . The hyperfine interaction is very 
weak in such systems. Due to the feedback, however, 
we determine a T c in the millikelvin range. The ordered 
phase leads furthermore to a universal reduction of the 
conductance and should be detectable by standard trans- 
port measurements. 



CONCLUSIONS 

We have discussed here the conditions necessary for fer- 
romagnetism of nuclear spins embedded in a 2DEG. 
Such a system is naturally in the RKKY regime and the 
interaction between the nuclear spins is carried by the 
electron spin susceptibility. Electron correlations are cru- 
cial, and the stability of the magnetic (ferro or different) 
phase depends on the nonanalytic behavior, linear in mo- 
mentum |q|, of the electron spin susceptibility. The non- 
analytic behavior cannot be found from standard Fermi 
liquid theory. It is a consequence of a strong renormaliza- 
tion of the electron-electron interaction which, to the best 
of our knowledge, strongly depends on a nonuniversal 
cutoff scale specific for the sample under investigation. 
The stronger the interaction, however, the higher also the 
critical temperature for the nuclear (ferro)magnet. Our 
estimates show that a transition temperature reaching up 
into the millikelvin range may be achievable. A simi- 
lar temperature range is estimated for one-dimensional 
conductors such as the recently available 13 C single-wall 
nanotubes, where the ordered phase should be detectable 
by conductance measurements. 

To conclude, we stress that such physics is not re- 
stricted to nuclear spins in metals. We expect a similar 
behavior for any Kondo-lattice system in low dimensions 
with interacting conduction electrons. 
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